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Abstract. Two-particle correlation functions of negatively charged hadrons from Pb- 
Pb collisions at 158 GeV/ c per nucleon have been measured by the WA97 experiment at 
the CERN SPS. A Coulomb correction procedure that assumes an expanding source 
has been implemented. Within the framework of an expanding thermalized source 
model the size and dynamical state of the collision fireball at freeze-out have been 
reconstructed as a function of the centrality of the collision. Less central collisions 
exhibit a different dynamics than central ones: both transverse and longitudinal 
expansion velocities are slower, the expansion duration is shorter and the system freezes 
out showing smaller dimensions and higher temperature. 
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1. Introduction 

The study of ultrarelativistic heavy ion collisions is motivated mainly by the QCD 
predictions that at sufficiently high energy density the excited nuclear matter undergoes 
a phase transition into a system of deconfined quarks and gluons (Quark-Gluon Plasma: 
QGP). 

The most direct measurement of space-time characteristics of the collision region 
is provided by Hanbury-Brown and Twiss (HBT) interferometry, a method developed 
originally in astronomy in the fifties |J and applied independently in the field of particle 
physics a few years later @. 

Over the last years this method has been widely used in ultrarelativistic heavy 
collisions to investigate the space-time geometry of the source, and ample evidence 
accumulated showing that information can be gained both on its geometry and its 
dynamics. For a recent review see e.g. ||. 

In ref. ||, Q an analysis strategy for reconstructing the freeze-out source state in 
heavy-ion collisions was suggested. A model for an emitting source characterized by 
a collective flow dynamics superimposed to a random (thermal) motion was assumed 
therein. The parameters of the model are the freeze-out temperature T, the transverse 
geometric (Gaussian) radius y/2R G , the average freeze-out proper time r , the mean 
proper time duration of particle emission At, the transverse flow velocity (3±, and Arj, 
that parametrizes the longitudinal extension of the source. 

This approach has been followed by the NA49 experiment, that has reconstructed the 
freeze-out state for very central Pb-Pb collisions (about 3% of the most central events) 
at 158 GeV per nucleon ||. 

In this paper we also follow the model of ref. ||, ||, but in addition we evaluate all 
the source parameters as a function of centrality, owing to the relatively large centrality 
window explored by the WA97 experiment. 

The correlation analysis presented here is based on pairs of negative hadrons 
(h~), which are dominated by pairs of identical pions. Misidentified particles lead 
to counting of unlike pairs which do not give rise to Bose-Einstein correlation. It has 
been shown in f?|, [8] (NA35 exp.) and in || (NA49 exp.) that the main effect of 
particle misidentification is to reduce the value of the chaoticity parameter A, which, 
however, is not used for the reconstruction of the size and dynamical state of the source; 
the relevant source parameters are effected by a few percents only. Moreover in this 
paper we present a compilation of recent results from other SPS experiments, which 
use either unidentified h'-h~ pairs (NA49 0, WA98 @) or Tr'-ir' pairs (NA44 fTO]); 
their agreement supports the assumption that HBT correlation can be performed using 
h~-h~ pairs. 

The paper is structured as follows: the WA97 apparatus is briefly described in 
section ^. The experimental derivation of the two-particle correlation function is 
discussed in section |||; section |] deals with the Coulomb corrections. In section [5] 
the measured radii are presented as a function of the pair momentum; from them 
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the source freeze-out size and temperature, the average freeze-out time and duration, 
and the expansion velocity are derived as functions of the collision centrality. From 
these observables a dynamical picture of the collision is depicted in section [|. Finally 
conclusions are drawn in section |7|. 



2. The WA97 experiment 

The WA97 set-up, shown schematically in fig. p], is described in detail in ref. The 
target and the silicon telescope were placed inside the homogeneous 1.8 T magnetic field 
of the CERN Omega magnet. The 158 A GeV/c lead beam from the CERN SPS was 



pad chambers 




Figure 1. The WA97 set-up. 



incident on a lead target with a thickness corresponding to 1% of an interaction length. 
Scintillator petal detectors behind the target provided an interaction trigger selecting 
approximately the most central 40% of the Pb-Pb collisions. Two planes of microstrip 
multiplicity detectors covering, for all pt values, the pseudorapidity region 2 < r\ < 3 
and 3 < rj < 4 respectively, provided information for more detailed off-line study of the 
centrality dependence of single and double particle spectra. In the present analysis the 
WA97 usual four centrality bins, described in ref. |L2], |13[ for the study of strangeness 
production, have been adopted. 

The core of the WA97 set-up is a silicon telescope consisting of 7 planes of pixel 
detectors [|14]] which have a pixel size of 75 x 500 /im 2 and of 10 planes of 50 fim pitch 
silicon microstrips. The telescope has 5 x 5 cm 2 cross section and contains ~ 0.5 x 10 6 
detecting elements. This precise tracking device was placed 60 cm behind the target 
slightly above the beam line and inclined by 48mrad, pointing to the target. 

The track reconstruction is done in the compact part of the telescope, where 11 
planes of silicon detectors are closely packed over a distance of 30 cm. The momentum 
resolution for fast tracks is improved using lever arm detectors, which consists of 
additional pixel and microstrip planes and of three MWPC's with a cathode pad readout. 
Due to the precise tracking and to its compactness the WA97 experiment is well suited 
for the measurements of the two-particle correlation functions. 
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3. Data analysis 

The quantum statistical correlation between two identical particles coming from a 
chaotic source yields information on the dimensions and the dynamics of the source. 
Experimentally the correlation is studied through a "correlation function" , defined as: 

<w«)=4g (i) 

where S(q) is the measured two-particle (Signal) distribution as a function of the relative 
four momentum q = p\—p 2 and B(q) is the reference (Background) distribution built by 
pairing two particles taken from different events. For each actual pair in one event (in the 
Signal) about 50 background pairs have been formed from events of similar multiplicity, 
so that the statistical error on C 2 is dominated by the signal. The normalization factor 
M is determined by imposing that the integral of 6*2(9) over a wide interval outside the 
region of sensitive correlation (which is at low q) is equal to one. 

From the correlation function C 2 the description of the source has been extracted 
using two different parametrizations. 

(1) The "standard" Cartesian (or Pratt-Bertsch) parametrization [[T5|, [L6| is expressed 
in the out-side-longitudinal (OSL) system; this has the longitudinal direction along 
the beam axis, while in the transverse plane the out direction is chosen parallel 
to the transverse component of the pair momentum K t , the remaining Cartesian 
component denoting the side direction. This parametrization is given by the 
following formula: 

C 2 (q) = 1 + A exp [-R 2 q 2 - R 2 s q 2 s - R 2 qf - 2\R ol \R oiqi q ] (2) 

where q D , q s , qi are the spatial components of q in the OSL system and R Q , R s , Ri, 
R i, which are fit parameters, provide information on the extension and dynamics 
of the source ||16|| . A is called the "chaoticity parameter" and is equal to unity for a 
completly chaotic source and smaller than unity for partially chaotic source. It is 
however affected in a non trivial way by resonance decays, particle misidentification 



and other effects [17], which are expected to depend on the pair momentum. 



(2) The Yano-Koonin-Podgoretskii (YKP) parametrization [18] uses the components 



Q± = yQo + 9s > <?o = Ei — E 2 , q\\ = qi and starts from the ansatz: 

C 2 (q) = 1+Aexp [-R 2 ± q 2 ± ~ 7jfc(?|| ~ v yk qo) 2 R\ ~ 7yk(lo ~ v yk q\\) 2 Rl] (3) 

where R±, R\\, Rq are called the "YKP radii", although R has temporal dimension 
(fm/c); v y k, called "Yano-Koonin velocity", is measured in units of c and 7^ = 
(1 " ^)- 1/2 - 

The transverse orientation of the reference system is different in the two 
parametrizations. Concerning the longitudinal direction, the LCMS (longitudinal co- 
moving system) frame, in which the average pair momentum has only a transverse 
component {K\ = 0), has been chosen for both of them. This system has the advantage 
of a more straightforward interpretation of the source parameters. Hence, a reference 
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Figure 2. Projections of the hr-hr correlation function measured using Cartesian 
(top) and YKP (bottom) parametrization for pair rapidities 2.55 < < 3.65 and 
pair trasverse momenta 0.23 < K t < 0.39GeV/c. Solid curves are the result of a 
least-square fit to a Gaussian function (see the text). 



system is used where the relative momentum q is measured on a pair-by-pair basis. The 
same is true in the transverse direction for the Cartesian parametrization. 



As claimed in ref. [I9[, fitting by the Least Square method and using the square 
root of the number of counts as an estimator of the error introduces a systematic bias. A 
Maximum Likelihood method has been applied by minimizing the negative logarithmic 
likelihood function: 

- 2\n£(R) =2j2[CiBi - ^hi(C^) + ln(^!)] (4) 

i 

where Cj is the correlation hypothesis for a given set of fit parameters R (e.g. R = 
(A, R±, R\\, Ro, v y k))- 

The main quality test of the full analysis procedure comes from the consistency 
between parameters obtained from Cartesian and YKP forms of the correlation function. 
This topic is discussed in section |[ Two secondary checks on the fit quality have been 
performed. Firstly, the x 2 /NDF-values are calculated for the parameters corresponding 
to the best likelihood fits and are found to be distributed around unity. Secondly 
we have fitted the projections of the 3-dimensional correlation functions onto one 
momentum difference component, with narrow cuts on the other components, using 
a least-squares method to a Gaussian function. This yields results consistent with the 
results of the 3-dimensional fits. A typical sample of these projections is shown in 
fig. § where the correlation function has been integrated along the other qj components 
in the intervals \qj\ < 20MeV, except for the q^ = Ei — E 2 projection where 
\q 3 \ < 50MeV/c(j=i_,||). 

The determination of the size and dynamical state of the source at freeze-out 
requires information about the dependence of the YKP radii on the mean momentum 
K = | (pi + P2) of the pair. This can be parametrized by its transverse component 
K t = W(Pyi+Py2) 2 + (p z i+p Z 2) 2 and the pair rapidity = \ log g±j2±jg±&a. 
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Figure 3. Two-particle acceptance window (assumed to be pions) used in the present 
analysis. The correlation function has been measured in each of the twenty rectangles 
patching the window together. The K t and average values, measured over the 
rectangles are reported in table 1. The dashed line is drawn at Y n7r = y cm . 



With this in mind, our (two-particle) acceptance window, which is shown in fig. 
has been binned in 20 rectangles in each of which the correlation function C^g) has 
been measured independently. Table 1 displays the average values of the K t and Y W7T 
distributions in each of the 20 bins used for the analysis. 



Table 1. Averages of the K t (uppermost number in the cells, in GeV/c units) and (bottom 
numbers) distributions in the 20 rectangles of fig. ||. 
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Additionally, in order to investigate the dynamics of the collisions as a function of 
centrality, the whole analysis has been repeated for each of the four centrality classes 
selected by WA97 to study the strange particle production JL2], [L3|]. 

The quality of a correlation measurement is determined by the resolution of the 
two-particle momentum difference ("relative momentum") q. The relative momentum 
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resolution was studied using a Monte Carlo chain based on the simulation code 
GEANT ||20|| . As reported in table 2, we have found that all relative momentum 
projections used in this analysis (q a , q s , q\\ qj_, q\\, go) are measured with an error not 
bigger than lOMeV/c over our two-particle acceptance window. 



Table 2. The resolution of relative momenta q ,q s ,qi and q±,q\\,qo in the LCMS frame. 
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The bin size used to build the correlation functions has been accordingly taken to 
be lOMeV/c in all relative momentum components. It has been shown in |2l| that 
due to the finite bin size the values of the extracted HBT radii are systematically 
underestimated, the greater the radius the higher the error. Since the typical half-width 
of the correlation function in our acceptance window is ~ 50 MeV/c (corresponding to a 
typical radius of about 4fm), our binning yields an underestimation smaller than 10%. 

The detector's ability to distinguish a pair of close tracks from a single track 
decreases with decreasing track separation, depending on the Si-pixel size. The 
correction for this inefficiency was calculated via a Monte Carlo simulation by 
reconstructing a sample of (generated) two-tracks events embedded in real ones. The 
maximum inefficiency was found to be ~ 4% for pairs populating the bin of the 
correlation functions centred at q = (relative to high q pairs). 

The results presented here are based on 17 million events collected in the 1995 run. 
Only tracks pointing to the main interaction vertex, which is reconstructed on a event 
by event basis, are selected. After event reconstruction we end up with 13 million h~ 
pairs in the correlation sensitive region \q\ < lOOMeV/c (LCMS system). 



4. Coulomb correction 



The observed two-particle correlation is the result of two different contributions, the 
Bose-Einstein effect and the Coulomb interactions. The Coulomb interaction between 
the particles of a pair accelerates them relative to each other, thus depleting the two- 
particle correlation function at small relative momenta. 

In order to get the Bose-Einstein correlation alone, the correlation function C% 
needs to be corrected for the Coulomb interaction among the particles, so that: 

C 2 (pi,P2) = C™ w { Pl ,p 2 ) X K C oul{pi,P2) (5) 
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where C™™ is the B-E correlation of the experimental data and Kc ou i is the correction 
for the Coulomb interaction. In the past, the correction was done by describing the 
emission of an isolated pairs of charged particles from a point-like neutral region (Gamow 
factor [E^, |23l]). This correction has been shown to be inadequate in heavy ion collision 



at the SPS energy [^4|, 2^, essentially due to the extent of the source which is far from 
a point-like approximation. The Gamow factor has been considered in this analysis just 
as a reference for the two correction methods implemented. Those are: 

i) The method based on the measure of oppositely charged particle correlation, fol- 
lowing the approach of Baym and Braun-Munzinger |2(| . 

Oppositely charged particles do not present correlation due to quantum statistical 
effects; they can thus provide a direct measure of the additional correlation intro- 
duced by the Coulomb interaction, accounting for the reverse in the charge sign. 
The method suggested in [26| is a classical model which neglects Coulomb inter- 
action inside a pair for separations less then a given r and includes it for larger 
separations. A charge screening effect, expected in the high multiplicity environ- 
ment created just after an ultra-relativistic heavy ion collision, is thus taken into 
account. The only free parameter of the model, vq, can be obtained by a fit to 
the measured data points of the correlator for oppositely charged particles. The 
obtained value, r = 6.1 ± 0.5 fm, is consistent with the NA49 measurement 0. 
The resulting correction function is shown in fig. f| as a function of qi nv = \J q 2 — 



ii) The full Coulomb wave-function integration This method has the 

advantage of being adaptable to any emission mechanism and can be implemented 
in bin-by-bin correction, independently of which particular projections of q = P\—p2 
are in use. Coulomb interactions with the residual nuclear system are neglected. 
As a validation test of the procedure, we first assumed a point-like source (i.e. very 
small volume, << lfm) with sudden emission (i.e. the particles are emitted at the 
same time), comparing to the Gamow correction. The agreement is satisfactory as 
shown in fig. [|. Then two kinds of source have been implemented: 

(a) a Gaussian source with sudden emission, which has been traditionally used in 
connection with the wave-function integration. 

(b) the same source model adopted to investigate the features of the expansion 
dynamics and geometry of the fireball. 

Since the parameters of the model, namely T, /3±, Rg, At/, At, r , are required 
to describe the source during the computation of the Coulomb correction, an 
iterative procedure has been followed: 

- a first complete analysis of the correlation functions corrected with the 
Baym and Braun-Munzinger approach, leads to a first estimate of the 
model parameters; 

- these parameters are used to calculate a new Coulomb correction by means 
of the full wave-function integration method; 
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- a second full analysis of the correlation with the new Coulomb correction 
(that will be presented in the next section) leads to a second, better 
evaluation of the model parameters. 
The procedure converged satisfactorily after the second iteration. 

A final comparison of all the methods is shown in fig. |j. As expected, the Gamow 
factor significantly overestimates the Coulomb interaction in heavy ion collision; the 
corrections calculated with the static Gaussian source and those of the Baym model, 
that implicitly also assumes a static source, are found to be very similar. The expanding 
source has proven to yield the smallest corrections. We may conclude that treating a 
rapidly expanding source as if it were a static system can lead to a slight overestimate of 
the Coulomb corrections. The correction has been applied on a bin-by-bin basis, acting 

















Gamow 
Baym (h"- 
Point like 
Gaussian 


h + ) 














• 




• 










1 


Expandin 


g 






■ 




















-1 . 





































l) 10 00 00 -10 00 60 70 30 90 0C 

q mv (MeV/c) 



Figure 4. Corrections to the correlation function for the Coulomb interaction 
according different models: the Gamow factor (dashed line), the Baym method based 
upon the correlation of oppositely charged particles (solid line) , the full wave- function 
integration for a point-like source (black squares) and for a (static) Gaussian source 
with rms radius of 4 fm (open circles) , the full wave-function integration for an 
expanding thermalized source ^ (shadow area). 

on the qi nv of that bin, even if we are making a correlation function in some other set 
of variables. 



5. Correlation study 

5.1. Cartesian and YKP parametrizations: consistency 

In the context of the model by Chapman et al. [|j] for a finite expanding thermalized 
source, the YKP parameters give more direct access to the features of the source. 
Therefore in this paper we shall follow this approach. However, the Cartesian 
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radii have been systematically used to test the goodness of the fitting procedure 
of the YKP correlation function, by means of equations relating the two sets of 
parameters (for these equations see e.g. [|, |29|). Fig. [| shows typical fit results for 
the parameters (A, R±, R», Rq) versus K t as obtained directly from YKP parametrization 
(black triangles) and deduced from the corresponding Cartesian set of parameters (open 
circles); as a second representative sample, in fig. [] are displayed the Cartesian radii 
versus Y nn evaluated directly and from the corresponding YKP correlation functions. It 
is worth stressing that the analysis has been carried out in a completely independent way 
for the two parametrizations: this means, for instance, that the corrections for Coulomb 
interaction and two-track resolution have been calculated and applied independently. In 
short, Cartesian and YKP radii give two separate measurements of the same physical 
quantities; their full compatibility gives confidence on the absence of systematic errors 
in the analysis. 
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Figure 5. The A factor and the YKP radii as calculated directly (black triangles) 
and from the Cartesian parametrization (open circles) as a function of K t in the full 
y w7r and centrality interval. The remaining parameter v y k is displayed in fig. ^ (plot 
labelled "ALL"). 



In fig. we have compiled recent results from other CERN heavy-ion experiments, 
namely NA49 |, ||, |0|, WA98 @ and NA44 0, along with our results for the two 
most central collision classes (III and IV) merged together. The overall agreement of 
our data with the other experiments is convincing both in the common K t region and 
when extrapolating to small K t . WA97 is the sole experiment that explores the high K t 
region and our data suggest a flattening of the R±, R\\, R s and Ri behaviour with K t 
at high values. 
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Figure 6. Cartesian parameters as calculated directly (open circles) and from the 
YKP parametrization (black triangles) as a function of Y wv in the interval 0.23 < 
K t < 0.76 GeV/c. 



5.2. The Yano-Koonin velocity versus Y^: longitudinal flow 

Fig. |8] shows the YK rapidity Yyk = \ log rz- 3 ^ + Y W7T , calculated from the fitted 
parameter v y k, as a function of the pair rapidity Y n% (both displayed in the laboratory 
system) separately for the four centrality classes and in the inclusive case (class "ALL"). 
The Vyk fit parameter is often determined with a larger relative error than R±, R\\ or A 
(the same is true for Rq); therefore four data sets (in the K t variable) have been merged 
in order to obtain small errors on v y k- 

As explained in ref. |J the dependence of the YK rapidity on the pair rapidity 
should directly measure the longitudinal expansion of the source and cleanly separate 
it from its transverse dynamics. In particular, for a non-expanding source Yyk 
would be independent of Y wn (horizontal lines in fig. ||), while for the case of a 
longitudinally expanding, boost-invariant source these variables would be completely 
correlated (diagonal lines in the same figure). It has been shown in ref. [29|| that, as a 
consequence of the reduced thermal smearing at large K t , the linear relation Yyk — Ymr 
is better satisfied the larger the transverse momentum of the particle pair. Hence this 
linear relation cannot be generally interpreted as evidence for a true boost-invariant 
longitudinal expansion. The dependence of the YK velocity on K t is expected however 
to be of minor importance. 
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Figure 7. Dependence of the HBT radii on the transverse momentum of the pair for 
Pb-Pb collisions at SPS energy. 



The results of fig. |8| thus would suggest that the source expands longitudinally in 
a nearly boost-invariant way for the three most central collision classes and with much 
less intensity for the most peripheral one. We observe, however, an unexpectedly strong 
dependence of v y k (hence of Y Y k) on the transverse pair momentum K t , as discussed in 
the next section. 

5.3. The Yano-Koonin velocity versus K t : v y k positive and close to c at high K t 

Fig. shows the K t dependence of v y k in the inclusive case of all centrality classes 
merged together (class "ALL"), and in the four centrality classes separately. In order to 
make a quantitative comparison of the behaviour of v y k in the different classes, we have 
parametrized linear function of K t : v y k = a + bK t . The slope b is found to be 

similar for all the four centrality classes with a value of about 1.5 (GeV/c 2 ) -1 , while the 
value of the intercept a is different for class I (a ~ —0.8 c) with respect to classes II, III 
and IV (a ~ —0.6 c). The parametrization works satisfactorily up to K t ~ 0.8GeV/c, 
but the linear approximation becomes unphysical {y y k/c > 1) for K t > £=2 ~ 1 GeV/c, 
therefore a saturation of v y k up to a value close to the speed of lighthas to be supposed. 
Positive v y k in the LCMS means that the diagonal lines of fig. |^ are crossed. Points 
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Figure 8. Y n7T dependence of the YK rapidity for the four centrality classes and for 
the inclusive case (ALL). Class I corresponds to the most peripheral collisions and 
class IV to the most central ones. Filled circles are data, open circles are data reflected 
about mid-rapidity. 



in fig. |8| are measured in the wide interval 0.23 < K t < 0.75GeV/c, the average K t 
over this interval being however ~ 0.45GeV/c. Since v y k crosses the zero value for 
K t ~ 0.4 GeV/c, looking back at fig. || we may claim that the experimental points are 
close to the 45°-degree line mainly because of their (mean) K t values. We can conclude 
that the slopes of the plots versus Y^J S {SYS being either LAB or CMS), 

often encountered in literature, are influenced strongly by K t . Therefore no conclusive 
predictions on the longitudinal expansion can be drawn by looking solely at the Y 7rn 
dependence of the Yano-Koonin velocity. It is a fact, however, that we measure smaller 
Vyk for centrality class I: the a intercept has a larger negative value and the zero crossing 
occurs at higher K t for this class. Since the mean K t of the pairs is found to be the same 
over the four centrality classes, therefore it can be deduced a less intense longitudinal 
expansion for class I. 
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Figure 9. i^T t dependence of the YK velocity in the four centrality classes and 
the inclusive (ALL) one. Black triangles are v y k's measured directly in the YKP 
parametrization, open circles are deduced from Cartesian radii. 



5.4- The R± radius: transverse flow and temperature 

While longitudinal "flow" is not necessarily a signature for nuclear collectivity, because 
it may be due simply to incomplete stopping of the two colliding nuclei, transverse 
collective expansion flow can only be generated by the build up of a locally isotropic 
pressure component. The transverse radius parameter R± is invariant under longitudinal 
boosts and thus is not affected at all by longitudinal expansion. It is expected to drop 
as a function of M t = a/ vn\ + Kf if the source expands in the transverse direction || . 

In ref. [§, in the context of a source model characterized by longitudinal and 
transverse expansion, the following expression for R± has been derived: 

i 

2 

(6) 

where Rq is equal to the transverse geometric (Gaussian) radius of the source times 
1 / y/2, T is the freeze-out temperature and f3± the slope of the (linear) transverse flow 
velocity profile: /3j_(r) = /3±j^. A bidimensional fit of function (|6]) to the experimental 

data has led to the determination of Rq and The fit results are summarized in 
table 3. 



R±(K t ,Y^) = R G 



i + M t ^cosh(F yi ,-y; 
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Table 3. The transverse geometric parameter Rq and the ratio (3j_/T, T being the freeze-out 
temperature. 



Centrality 


Rg (fm) 


§- (GeV- 1 ) 


X 2 /ndf 


ALL 


4.75 ±0.38 


1.62 ±0.60 


21.6/18 


I 


3.20 ±0.31 


0.68 ±0.55 


21.6/12 


II 


4.63 ± 0.44 


1.82 ±0.78 


15.5/16 


III 


5.03 ±0.56 


1.97 ±0.92 


20.6/17 


IV 


5.06 ±0.58 


1.76 ±0.78 


23.0/16 



The K dependence of R± determines only the ratio ^r, so that a full range for the 
(/3±, T) doublet is still conceivable. WA97 has measured the inverse m t slope parameters 
T e ff as a function of centrality ||31|| . This slope can be thought as the inverse of a blue- 
shifted temperature T e jf = T-^ ^^ ". This relation provides a second constraint on 
the (/3j_,T) doublet range. Combining the two-particle correlation measurements with 
the single particle transverse mass spectra, disentangles temperature T and transverse 
flow P±. The allowed regions in the (T,{3±) parameter space, as determined by the 
combination of these two measurements, are shown in fig. [Tl], and the disentangled 
values for T and /3± are reported in table 4. 



Table 4. The freeze-out temperature T and the intensity of the transverse flow j3±. 



Centrality 



T(MeV) 



01 



ALL 
I 

II 
III 

IV 



+12 
-10 
+26 
-13 

121±lf 
11711? 

-15 
-11 



123 
140 



120 



n 44+0-06 

o so +0 - 09 

47+0- 07 

u -^' -0.10 

0.481°;?? 

h0.07 
-0.10 



0.46; 



Postponing to the next section a global discussion about these results, we wish 
here to emphasize that class I features a smaller transverse freeze-out radius, a lower 
transverse expansion velocity and a higher temperature. This suggests an expansion 
on a smaller scale; the higher temperature at the freeze-out may be interpreted as the 
remnant of an earlier decoupling of the expanding system. 

5.5. The R\\ radius: average freeze-out proper time 

In the same model of ref. ||], the dependence of Rn on (K t , Y nn ) determines the average 
freeze-out proper time To via the formula: 

R\\ = r [^ cosh(V yi , - - cosh- 2 (V y ^ - Y^) + l/^) 2 }^ (7) 
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Figure 10. Allowed regions in the freeze-out temperature vs. transverse flow velocity 
space (see text). Bands are drawn at ±a around fitted values. 



The parameter An is fixed by the extension of the single particle rapidity distribution || 



it has been measured by the large acceptance NA49 experiment to be An = 1.3±0.1 |32 
The results of the bidimensional fits of function (|7|) to the experimental data are 
presented in table 5. 



Table 5. Average freeze-out proper time tq of the source. 



Centrality 


r (fm/c) 


X 2 /ndf 


ALL 


5.35 ±0.20 


8.6/19 


I 


3.66 ±0.22 


24.3/13 


II 


5.08 ±0.23 


16.5/17 


III 


5.56 ±0.23 


20.9/19 


IV 


5.58 ±0.24 


9.6/17 



The centrality classes III and IV show a similar behaviour, being characterized by a 
common mean freeze-out (proper) time To — 5.6 fm/c, class II presents a sligthly shorter 
To of about 5 fm/c. All these values are much larger than the 3.7 fm/c measured for class 
I. 

5. 6. The -Rq parameter: mean duration of the emission 

The YKP parameter R Q is related to the mean emission duration At of the source |33[] : 
At ~ R at high K t (8) 
% Arj is related to the width of the rapidity distribution Ay by (Ay) 2 = (Arj) 2 + 
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As one can see for the inclusive case of ALL centralities shown in fig. [II], the Rq radius 
does not present a significant dependence on the pair rapidity Y n7T . Therefore we have 
merged all the four bins with the aim of reducing the statistical fluctuations on the 
extracted R when the four centrality classes are investigated separately. The Rq fit 



0.23SK <0.39 GeV/c 



0.39SK <0.49 GeV/c 



I 

rj 

of 



0.49<K,<0.59 GeV/c 



J 
o 



02_ 



0.76<K,<1.3 GeV/c 



0.59<K,<0.76 GeV/c 



Figure 11. The parameter as a function of Y n7T in five successive Kt intervals. 
The dotted horizontal lines through the points are the best-fit results with a constant 
function. The showed data refer to the case of "ALL" centralities. 



results for these data sets are displayed in fig. [L2| as a function of K t . 

Negative values of the Rq parameter can be expected from an opaque source even at 
K± = p4|. In an opaque source the freeze-out occurs in the form of surface emission, 



as in the case of photons radiated by the sun. Fig. [L2| suggests that all centrality classes 
are compatible with a transparent source that suddenly freezes out (At ~ 0). 
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Figure 12. Dependence of -Rq on K t in our full rapidity interval, for the different 
classes of collision centrality. 
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The study of the i^-dependence of YKP parameters, combined with single-particle 
transverse-mass spectra, has determined several features of the source: its freeze-out 
temperature T, its transverse geometric (Gaussian) radius \f2Rc-, its average freeze-out 
proper time r , the mean proper time duration of particle emission At, the strength 
of the transverse flow /3±, some hints on the longitudinal expansion. In table 6 we 
summarize the main results of this study. These data provide implicitly a dynamical 
picture of the collision process. 

Table 6. Source parameters for the four centrality classes. 





IV 


III 


II 


I 


Rg (fin) 


5.1 ±0.6 


5.0 ±0.6 


4.6 ±0.4 


3.2 ±0.3 


T(MeV) 


i2oi j; 


ii7+;« 


121+15 


1401^ 


P±(c) 


46 +U U7 




47+"- 07 

-0.10 


30+°' 09 


T (fm/c) 


5.6 ±0.2 


5.6 ±0.2 


5.1 ±0.2 


3.7 ±0.2 


Ar(fm/c) 


w at high K t 


Long. 




fast 




slow 


expans. 











With the aim of seeing if this picture is self-consistent, we can compare the two- 
dimensional rms width R(.^ zeout = \/2Rq with the two-dimensional rms width of a cold 
lead nucleus 

iC,= y|l.2A*^4.5fm. (9) 

For most central collisions (class III and IV), the system expands by a factor ~ 1.6 
or, equivalently, by 2.65 fm in the transverse direction. If the transverse flow velocity 
is equal to (3± ~ 0.47 c during the whole expansion, in a time of tq ~ 5.6fm/c nuclear 
matter would travel over Tof3± ~ 2.6 fm in the transverse direction. This is consistent 
with the previous estimate from the difference R{.^ zeout — R^ms- Therefore the HBT 
description is internally consistent. 

For less central collisions (classes I, II) the initial overlap of colliding nuclei is 
expected to be smaller than that given by eq.( [|). A realistic estimation of such an 
overlap is not straightforward. However, a backward reconstruction from the data is 
still feasible. If the freeze-out transverse two-dimensional rms width is 4.5 fm (6.5 fm) 
for class I (class II), accounting for an overall expansion of Tof3± = 1.1 fm (2.4 fm), we 
may estimate the initial transverse rms width to be 3.4 fm (4.1 fm). 

No quantitative description of the longitudinal expansion can be done. We can 
qualitatively say that the longitudinal expansion becomes stronger the more central the 
collisions; additionally our data would suggest a sudden step passing from centrality 
class I to class II (see fig. 

Concerning the collision time development, we have argued that the pion emission 
is a sudden process (At ~ 0), starting after ~ 5.5fm/c (3.7fm/c) in the case of most 
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central (peripheral) collisions, that happens in bulk (i.e. not surface dominated). This 
resembles the decoupling process of photons in the Early Universe. 

It can be supposed that pion decoupling is caused mostly by a rapid cooling 
and dilution of the baryon density. In thermal models |35j the chemical freeze- 
out temperature (~ 170 MeV), needed to describe the observed particle ratios, is 
considerably higher than the thermal freeze-out temperature. This means that the 
relative abundances of particle species stop changing earlier than the end of (strong) 
elastic collisions (thermal freeze-out). The last stage of expansion would thus be 
characterized by a volume dilution of baryon density without any appreciable particle 
recombination. 

7. Conclusions 

The study of two-particle correlation functions has provided valuable information both 
on the geometry and the dynamical state of the reaction zone at freeze-out. A Coulomb 
correction procedure that assumes the same model for pion emission as that used in the 
HBT analysis has been implemented. 

We have measured both negative and positive Yano-Koonin velocities in the LCMS 
system at slightly forward rapidity (with respect to y cm ), observing a strong increase of 
Vyj, with the transverse pair momentum K t . 

The simultaneous analysis of single-particle transverse mass spectra and two- 
particle correlation, with the help of a model that accounts for both longitudinal and 
transverse expansion of the source, has allowed for a reconstruction of the final freeze- 
out state of the source. We have measured the following parameters as a function of 
centrality: 

• the thermal freeze-out temperature, 

• the rms transverse Gaussian radius, 

• the average freeze-out proper time, 

• the mean proper time duration of particle emission, 

• the transverse flow velocity. 

From these data we have shown that our less central collision class (I) exhibits a different 
dynamics: the indications are that the whole expansion develops less intensely (both 
longitudinally and in transverse direction), the duration of expansion is shorter, the 
final (i.e. at freeze-out) source dimension is smaller, the system ends up with a higher 
temperature than in the case of more central collisions. The source parameters of the 
two most central collision classes (III and IV) have been found to be very similar, 
thus suggesting the same geometrical and dynamical freeze-out state for these classes. 
Class II presents parameters that are intermediate between class I and class III-IV, its 
dynamics being however much closer to that of classes III-IV. 
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